Study of the de Almeida-Thouless line using one-dimensional power-law diluted 

Heisenberg Spin Glasses 
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We test for the presence or absence of the de Almeida-Thouless line using one-dimensional power- 
law diluted Heisenberg spin glass model, in which the rms strength of the interactions decays with 
distance, r as l/r"^. It is argued that varying the power a is analogous to varying the space 
dimension d in a short-range model. For a = 0.6, which is in the mean field regime regime, we find 
clear evidence for an AT line. For a — 0.85, which is in the non-mean-field regime and corresponds 
to a space dimension of close to 3, we find no AT line, though we cannot rule one out for very 
small fields. Finally for a = 0.75, which is in the non-mean-field regime but closer to the mean-field 
boundary, the evidence suggests that there is an AT line, though the possibility that even larger 
sizes are needed to see the asymptotic behavior can not be ruled out. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 



I. INTRODUCTION 



One of the major unsolved questions in the field of spin 
glasses is the nature of the low-temperature glassy phase. 
Two main candidate theories exist: the replica symme- 
try breaking (RSB) theory which assumes that real spin 
glasses behave in a similar way to the mean-field theorjti, 
and the droplet model^ii^, which provides a phenomeno- 
logical approach. An important difference between the 
predictions of the two approaches is whether or not there 
is a line of phase transitions, known as the de Almeida- 
Thoules8i (AT) line, in the presence of a magnetic field. 
An AT line is predicted in the RSB picture, but is argued 
not to exist according to the droplet picture. The exis- 
tence of the AT line is arguably the most striking features 
of the immensely complicated mean-field-theory of spin 
giassesi, and so it is of intrinsic interest to know whether 
or not it occurs in real spin glasses. Some time ago Bray 
and Roberta^ investigated whether or not there is an AT 
line below six dimensions (six is the upper critical dimen- 
sion for spin glasses) using renormalization group ideas, 
but did not find a stable, accessible fixed point. More 
recently, Temesvari^ claimed to be able to follow the AT 
line down to just below six dimensions. 

Theoretically, it is interesting to understand spin 
glasses in a range of dimensions. However, this is hard 
since one cannot simulate systems with a large number 
of spins N, where N = L''', for a range of linear sizes 
L (needed to do finite-size scaling) if the dimension d is 
large. In particular, the the mean-field regime, d > 6, 
is not directly amenable to simulation. Hence, instead, 
we investigate a one- dimensional model in which the rms 
strength of the interactions fall off with a power a of the 
distance. This model has the advantage of allowing one 
to study large (linear) sizes. Furthermore, it has been 
suggested that changing the value of a is analogous to 
changing the value of the space dimension d in a short- 
range model. Consequently we can study models in both 
the mean-field and non-mean-field regimes. The nature 



of the spin glass phase in long-range models was discussed 
recently by Moored. 

Most previous work^rJ^ on one-dimensional, long- 
range spin glass models, used Ising spins. However, in 
this paper, we build on a recent work by us^^ which 
showed that, in mean field theory, an AT line also occurs 
in m-component vector spins provided the magnetic field 
is random in direction. More precisely, we perform Monte 
Carlo simulations on the three-component (Heisenberg) 
spin glass in one dimension with long-range interactions 
in the presence of a random magnetic field. This work 
follows on from a recent paper— where we considered the 
same model in zero-field. 

The plan of this paper is as follows: in Sec. [IT] we de- 
scribe the model and the Monte Carlo method used to 
simulate it, in Sec. lHII we describe the results, and finally 
in Sec. II VI we summarize our conclusions. 



II. MODEL AND METHOD 



The Hamiltonian we study is 



H ^ — Jij^i ■ Sj — hi ■ Si 



(1) 



where the Si, i — 1,2, ■ ■ ■ , N, are classical 3-component 
Heisenberg spins of unit length, and the interactions Jij 
are independent random variables with zero mean and a 
variance which falls off with a power of the distance Vij 
between the spins. 
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(2) 



The notation [• • • ]av indicates an average over the 
quenched disorder. In addition we set Ja — 0. The 
magnetic fields /if, where fj, denotes a cartesian spin com- 
ponent, are chosen to be independent Gaussian random 
fields, uncorrelated between sites, with zero mean, which 
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satisfy 



We set J = 1 so that, for the situation here, 



[/if /if]av = hi ^ij ^tJ 



(3) 



J ~ 1, z = 6, m ~ 3 , 



(9) 



Following Leuzzi et. al.— , and continuing along the lines 
of the zero- field paper—, the interactions of our model are 
such that, instead of the magnitude of the interaction 
falling off with distance like Eq. ([2]), it is the probability 
of there being a non-zero interaction between sites 
which falls off, and when an interaction does occur, its 
variance is independent of . The mean number of non- 
zero interactions from a site, which we call z, can be 
fixed, and here we take z — 6. To generate the set of pairs 
(z, j) that have an interaction with the desired probability 
we choose spin i randomly, and then choose j i) at 
distance rij with probability 



-2(7 



Ptj = 



(4) 



where, for r^j, we put the sites on a circle and use the 
distance of the chord, i.e. 



N 



(5) 



If i and j are already connected, we repeat the process 
until we find a pair which has not been connected before. 
We then connect i and j with an interaction picked from 
a Gaussian interaction whose mean is zero and whose 
standard deviation is J, which we set equal to 1. This 
process is repeated precisely Ni, — zN/2 times. 

The result is that each pair will be connected 

with a probability Pij which must satisfy the condition 
Pij = since Pij only depends on |i — j|, Pa ~ 0, 
and there are precisely Nz/2 connected pairs. It follows 
that, for a fixed site i, 



[ Jij ]av — J Pij — J 



(6) 



The mean-field spin glass transition temperature for m- 
component vector spins is given, for zero field, hy^ 



nMF 



( 7 , [ •Jij ]av I 
j J 



1/2 



■ J, 



(7) 



where the last equality follows from Eq. ([5]). In 
mean-field theory, the critical magnetic field at zero- 
temperature for m-component vector spinsH is (after ac- 
counting for the different normalization of the spins in 
Ref. [131) given by 



MF 



m{m — 2) 



(8) 



we have 



T, 



MF 



0.816. 



(10) 



the same as for the nearest-neighbor Heisenberg spin 
glass on a simple cubic lattice, and 



V2~ 1.414. 



(11) 



The ratio of these two quantities, which we shall refer to 
later, is given by 



h: 



MF 



j^MF 



V3~ 1.732. 



(12) 



According to Ref. [13|, a good approximation for the 
AT line in the mean field theory of the Heisenberg (m = 
3) spin glass is 



uMF 
"-AT 



(T) 



rpMF 
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m + 2 
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where 
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In Eq. ([T5)) we have again allowed for the different nor- 
malization of the spins in Ref. [13|. Equation ((T5)) is ex- 
act, in mean field theory, near T*^^, and, for m = 3, 
works very well down to quite low temperatures, see 
Fig. 1 in Ref. [3. Even for T = 0,(< = 1), Eq. ^ 
gives hf^/T^^^ = y/12/5 ~ 1.549, whereas the cor- 
rect value in mean field theory is, according to Eq. (|12l) . 
hfF /T^iP 1.732. 

We perform Monte Carlo simulations for this model 
in a magnetic field for three values of u: 0.6, 0.75, and 
0.85. As discussed in the zero-field paper—, 0.6 lies in 
the mean-field regime, while the other two values of a 
are in the non-mean-field regime. 

An approximate connection between a value of a 
and the effective dimension of an equivalent short-range 
model, doff, 



(15) 



2a- 1 



In the non mean-field region, a more accurate connec- 
tion, which involves the exponent 775^ of the short-range 
model, can also be obtained^^, but we will neglect this 
correction here since we have little information on rjsR. 
From Eq. (jlSp the effective dimensions corresponding to 
(T = 0.75 and 0.85 are dcs = 4 and d^s ~ 3, respectively. 

We continue to use the technology described in 
the zero-field paper—: overrelaxation sweeps, heatbath 
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sweeps, and parallel tempering. We perform one heat- 
bath sweep, and one parallel tempering sweep for every 
ten overrelaxation sweeps. 

The Gaussian nature of the interactions and the mag- 
netic fields affords a useful test for equilibration. The 
relation 



U : 



J2 Z , 



qs) + ^ (<z ■ 



1), 



(16) 



is valid in equilibrium but, very plausibly, the two sides 
approach their common equilibrium value from opposite 
directions as equilibrium is approached. Here 



U = - 



(17) 



is the average energy per spin, q = {i/N)J2i[{^i) ' 
(Si)]av is the Edwards-Anderson order parameter, qi = 
Si ■ Sj)^]av is the "link overlap", and 
:((S, • S,)2)]av, where = zN/2, 



and Cij = 1 if there is a bond between i and j and is zero 
otherwise. Equation (jl6p is easily derived by integrating 
by parts Eq. (|17p with respect to Jij and since they 
have Gaussian distributions. 

We determine both sides of Eq. p6)) for different num- 
bers of Monte Carlo sweeps (MCS) which increase in a 
logarithmic manner, each value being twice the previ- 
ous one. In all cases we average over the last half of the 
sweeps. We consider the data to be equilibrated, if, when 
averaging over a large number of samples, Eq. (|16|) is sat- 
isfied for at least the last two points. Note that in the 
numerics we set J = 1 . Table U lists the parameters of 
the simulation. 

We determine the wave-vector-dependent spin-glass 
susceptibility, given hyiS, 



Xscik) 



where 



N ^ m 



Jk (i-j) 



isms';), 



(18a) 



(18b) 



in which (■ ■ ■ ) denotes a thermal average and [• • ■ ]av an 
average over the disorder. To avoid bias, each thermal 
average is obtained from a separate copy of the spins, so 
we simulate four copies at each temperature. The spin 
glass correlation length is then determined from 



1 



f Xsg(O) 

2sin(fc,„in/2) V 

XSG (^inin) 



1 



1/(2<T-1) 



(19) 



where /cmin = {2tt/N). 

According to finite-size scalingi^iii, the correlation 
length of the finite-system varies, near the transition tem- 
perature Tc, as 



TV/a 



1^X[N^/-{T~T,)] , (2/3<(7< 1), 

jl/3 



= X[N^^-\T - T,)] , (1/2 <a< 2/3), 



(20a) 
(20b) 



TABLE I; Parameters of the simulations. A'^samp is the num- 
ber of samples, NcquH is the number of overrelaxation Monte 
Carlo sweeps for equilibration for each of the 2Nt replicas 
for a single sample. The same number of sweeps is done in 
the measurement phase, with a measurement performed every 
four overrelaxation sweeps. The number of heatbath sweeps 
is equal to 10% of the number of overrelaxation sweeps. Tmin 
and Tmax are the lowest and highest temperatures simulated, 
and Nt is the number of temperatures used in the parallel 
tempering. 



a 


hr 


iV 


samp 






Tin ax 


Nt 


0.6 


0.1 


128 


6000 


512 


0.20 


0.50 


21 


0.6 


0.1 


256 


6000 


512 


0.20 


0.50 


23 


0.6 


0.1 


512 


6000 


2048 


0.20 


0.50 


24 


0.6 


0.1 


1024 


6000 


4096 


0.20 


0.50 


27 


0.6 


0.1 


2048 


6000 


8192 


0.20 


0.50 


30 


0.6 


0.1 


4096 


2000 


16384 


0.25 


0.45 


40 


0.6 


0.1 


8192 


1000 


32768 


0.27 


0.45 


42 


0.6 


0.1 


16384 


1000 


16384 


0.30 


0.42 


45 


0.75 


0.1 


128 


6000 


8192 


0.05 


0.24 


21 


0.75 


0.1 


256 


6000 


16384 


0.05 


0.24 


22 


0.75 


0.1 


512 


6000 


16384 


0.10 


0.30 


24 


0.75 


0.1 


1024 


3500 


65536 


0.10 


0.24 


20 


0.75 


0.1 


2048 


2300 


131072 


0.10 


0.24 


23 


0.75 


0.1 


4096 


1000 


262144 


0.10 


0.30 


45 


0.85 


0.05 


128 


10000 


32768 


0.015 


0.20 


25 


0.85 


0.05 


256 


9000 


65536 


0.015 


0.20 


34 


0.85 


0.05 


512 


5000 


262144 


0.02 


0.20 


38 


0.85 


0.05 


1024 


5000 


524288 


0.03 


0.20 


40 


0.85 


0.05 


2048 


4100 


524288 


0.05 


0.20 


34 


0.85 


0.1 


128 


6000 


8192 


0.05 


0.25 


21 


0.85 


0.1 


256 


2500 


16384 


0.05 


0.25 


22 


0.85 


0.1 


512 


2000 


32768 


0.05 


0.25 


24 


0.85 


0.1 


1024 


1200 


262144 


0.05 


0.20 


34 


0.85 


0.1 


2048 


1600 


262144 


0.05 


0.20 


40 


0.85 


0.1 


4096 


1800 


131072 


0.10 


0.20 


30 



in which v, the correlation length exponent, is given, in 
the mean-field regime, by = l/(2cr— 1). It follows that, 
if there is a transition at T = T^, data for ^/N {S^/N"/^ in 
the mean-field region) for different system sizes A'' should 
cross at Tr. 



We also present data for xsG 
finite-size scaling form 



Xsg(O), which has the 



XsG 

XSG 
iVl/3 



= C[N^'''{T - TJ] , (2/3 <<7<l), (21a) 
= C[N^'^{T - TJ] , (1/2 < cr < 2/3). (21b) 



Hence curves of xsg/A'^"'' (xsg/A'"^/^ in the mean-field 
regime) should also intersect. This is particularly use- 
ful for long-range models since rj is given by the simple 
expression 2 — rj = 2a — 1 exactly. 

In practice, there are corrections to this finite-size- 
scaling, so data for different sizes do not all intersect 
at the exactly the same temperature. Including lead- 
ing corrections to scaling, the intersection temperature 
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0.01 




FIG. 1: (Color online) A finite size scaling plot of data for ^SG, 
for a = 0.6 in a magnetic field of hr = 0.1. The intersections 
indicate that there is a spin glass phase transition (AT line) 
in the presence of a magnetic field. 



FIG. 2: (Color online) A finite size scaling figure of data for 
XSG, for a = 0.6 with hr = 0.1. As in Fig.[Tl the intersections 
indicate the presence of an AT line. 



0.42 



T*{N, 2N) for sizes N and 2N varies a^^-^^ 

A 



T*{N,2N) =Tc + 



(22) 



where A is the amplitude of the leading correction, and, 
in the non mean-field regime, the exponent A is given by 



A = i + w 



(23) 



where uj is the leading correction to scaling exponent. 



III. RESULTS AND ANALYSIS 

A. cr = 0.6 

We recall that a = 0.6 lies in the mean-field regime. 
In this regime, simulations of the corresponding Ising 
model^ii^ found an AT line. From our plots for the 
Heisenberg spin glass in Figs. [T] and [5] for hr = 0.1, 
we come to the same conclusion here. According to 
Eqs. (I20bl) and (l2lbl) . data for S.sg/N"/^ and xsg/N^^^ 
should intersect at the transition temperature. We do, 
indeed find intersections, though the intersection temper- 
atures vary somewhat with size. 

The intersection temperatures are shown in Fig.|31 Fit- 
ting the intersection temperatures to Eq. ((22|) , using the 
known value^.^ A = 0.467 we find Tsg = 0.406 ± 0.012 
from S^sG (omitting the two smallest sizes), and Tsg — 
0.405 ± 0.007 from xsG (including all the data). These 
two results agree well with each other. Note that the 
intersection temperatures increase with increasing size, 
which suggests that they will not disappear in the ther- 
modynamic limit. 
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FIG. 3: (Color online) A plot of the intersection temperatures 
T'(N,2N) for a = 0.6, obtained from the data in Figs. [T] 
and[2l as a function of N^'^, with A fixed to its value 0.467, 
which is known exactly because a — 0.6 is in the mean-field 
regime. The fits give: Tsa = 0.406 ± 0.012 from ^sG and 
Tsg ~ 0.405 ±0.007 from xsg- For comparison, the zero field 
transition temperature ia^ Tsg — 0.563. 



It is interesting to compare this point on the AT line, 
{T,hAT{T)) = (0.405,0.1), with mean field predictions. 
Replacing T^^^ with the actual zero field transition tem- 
perature oSl4 Tc = 0.563 in Eqs. ([13]) and dH]), we find 
that T = 0.405 gives a field hAT{T = 0.405) = 0.130, 
which is slightly larger than the actual field value of 0.1. 
Hence the value of the field on the AT line for a = 0.6 is 
somewhat less than that expected in mean field theory, 
even allowing for the reduction in Tc from its mean field 
value. 

To conclude this section, the data suggests that there 



5 




0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 ..-X 



FIG. 4: (Color online) A finite size scaling plot of data for 
^SG for a — 0.75. The magnetic field is hr — 0.1. 




0.1 I 1 1 1 1 1 '-^ 1 

0.1 0.12 0.14 0.16 0.18 0.2 0.22 0.24 

T 

FIG. 5: (Color online) A finite size scaling plot of data for 
XSG for a — 0.75. The magnetic field is hr — 0.1. 



is an AT line for a — 0.6, as found earlier^'^^ for the 
corresponding Ising model. 

B. a = 0.75 

According to Eq. p3|) (which is approximate) this 
value of a corresponds to a short-range model in d = 4 
dimensions. For the corresponding Ising study, Ref. (l2j 
did not find an AT line, though this conclusion was sub- 
sequently challenged in Ref. 9] . 

Data for the spin glass correlation length and suscep- 
tibility are shown in Figs. |3] and [5] There are clearly 
intersections. Since we are now in the non mean- field 
regime, the exponent A in Eq. (12^ is not known and so 
should be treated as a fit parameter. However, the in- 
tersection temperatures have quite large error bars, and 



FIG. 6: (Color online) A plot of the intersection temperatures 
r*(Af, 2N) for a = 0.75, obtained from the data in Figs.[4]and 
O as a function of A'^"^ with A — 0.44, the value found for 
a = 0.75 in our zero field stud}»i^. We were not able to deter- 
mine A with any degree of precision from the data in Figs. |4] 
and [5] Both sets of data are compatible with a transition 
temperature of about 0.16 though there is also a hint in the 
data that the intersection temperatures start to drop for the 
largest sizes. For comparison, the zero field transition tem- 
perature isi^ Tc ~ 0.357. 



do not seem to vary monotonically, as shown in Fig. |6l 
Thus the data is not of good enough quality to determine 
A, and, in the plot, we have, rather arbitrarily used the 
value A = 0.44 obtained for a — 0.75 in our zero-field 
studyi\ 

Both sets of data are compatible with a non-zero Tc of 
about 0.16, which is to be compared with the zero field 
spin glass transition temperature Tsg — 0.357. However, 
there is a suggestion in the data that the values of T* 
decrease for the largest pairs of sizes, so perhaps one 
should treat this conclusion with some caution. 

As we did above for a = 0.6, we compare the pu- 
tative point on the AT line, {T,hAT{T)) = (0.16,0.1), 
with mean field predictions. Replacing with the 

actual zero field transition temperature of^ Tc = 0.357 
in Eqs. ^ and HH), we find that T = 0.16 gives a field 
HatIT — 0.16) — 0.227, which is considerably larger than 
the actual field value of 0.1. Hence, if the intersections in 
Figs.[l]and[51 do represent a transition in a field, the value 
of this AT field is is considerably less than that expected 
in mean field theory, even allowing for the reduction in 
Tc from its mean field value. 



C. cr = 0.85 

According to Eq. p^ . a — 0.85 corresponds to a short- 
range model in close to three dimensions. In their study 
of the Ising version of this model, Ref. [l3 did not find 
an AT line for this value of a. Ref. 0] did not consider 
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FIG. 7: (Color online) A finite size scaling plot of data for 
CsG for a = 0.85. The magnetic field is hr = 0.1. No in- 
tersections are found for the range of temperatures that we 
simulate indicating that there is no phase transition in this 
range. For comparison, the zero field transition temperature 
igii ~ 0.165. 



FIG. 9; (Color online) A finite size scaling plot of data for 
^SG for a — 0.85. The magnetic field is hr — 0.05. The data 
appears to merge at the lowest temperatures. The inset shows 
an enlargement of a region of the data. 
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FIG. 8: (Color online) A finite size scaling plot of data for XSG 
for a = 0.85. The magnetic field is K = 0.1. At T ~ 0.05, 
the largest sizes for which we have data. A'' = 512, 1024 and 
2048 merge together. For comparison, the zero field transition 
temperature ia^ Tc — 0.165. 



this value of a. 

We show finite-size scaling plots for S^sG and xsG for 
hr = 0.1 and 0.05 in Figs. EHTOl For hr = 0.1, the 
^SG data does not intersect down to the lowest temper- 
ature, 0.05, which is to be compared with the zero field 
transition temperaturei^ Tc ~ 0.165. The xsG data for 
hr = 0.1 seem to merge at the lowest temperature for the 
largest sizes that could be studied in this region. 

For hr = 0.05, the data for both £^sG and xsG merge 
at the lowest temperatures. One possible explanation 
of this is that the critical field at T = is /ic ^ 0.05. 
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FIG. 10: (Color online) A log-log finite size scaling plot of 
data for xsG for a — 0.85. The magnetic field is hr = 0.05. 
The data appears to merge at the lowest temperatures. 



Another possibility is that, for this small field, we are 
in a crossover region between zero-field behavior, where 
there is a transition, to finite-field behavior where there 
is none. 

According to Eq. in mean field theory the zero- 

temperature critical field, he is 1.732 times the zero- 
field transition temperature Tc- If we assume, based on 
our data, that he ^ 0.05 then, using Tc ~ 0.165 from 
Ref. [l4|, we have a ratio hc/Tc ^ 0.30, which is about 
17% of the mean field result, i.e. considerably smaller. 
Hence, if the data in Figs. |9] and [10] is interpreted to 



show a critical field of around 0.05, then this is much 
smaller than in mean field theory, even allowing for the 
(substantial) reduction in the zero-field Tc relative to the 
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mean field prediction. 

IV. CONCLUSION 

We have studied existence or otherwise of the AT line 
in the 1-dimensional Heisenberg spin glass with interac- 
tions which fall off as a power of the distance. We are 
able to study a large range of sizes in the temperature 
range of interest: N < 16384 for a = 0.6, N < 4096 for 
a = 0.75, and N < 2048 for a = 0.85 (up to 4096 at 
somewhat higher temperatures). 

For a = 0.6, which is in the mean- field regime (cr < 
2/3), we find an AT line. For a — 0.75, which is in the 
non-mean-field regime and corresponds to a short-range 
model with dimension about 4, the data does appear to 
find a phase transition for a field hr = 0.1, though the 
intersection temperatures in the finite-size scaling plots 
drop for the largest sizes, which might give one pause 
to accept that conclusion with certainty. For a — 0.85, 
which corresponds to a dimension of about 3, the data 
for hr — 0.05 can be interpreted as indicating that this 
value of field is close to a critical AT field in the limit 
of zero temperature. However, it can also be interpreted 
as indicating a crossover between the zero-field transition 
and behavior in a field which has no transition. 

For the corresponding Ising model, Ref. [l^l finds an 



AT line in the mean-field regime but not in the non-mean- 
field regime. This conclusion was challenged in Ref. Q, 
who claim that there is also an AT line in the non-mean- 
field regime, at least up for a between 2/3 and 0.75. A 
motivation for the present study was to see if a clearer 
numerical picture can emerge from the Heisenberg spin 
glass, where it is possible to study larger sizes than for 
the Ising case. Unfortunately, it seems that corrections 
to scaling are quite large in the Heisenberg case, so we 
are not able to give a precise value for the lower critical 
dimension of the AT line from the data in our paper. If 
there is no AT line, the system breaks up into domains 
of size £ (Imry-Ma length) which can be large at low 
temperatures. A possible explanation of our results is 
therefore that e{T 0) for a = 0.75, K = 0.1 is larger 
than the largest system size, N — 4096, and that for 
(T = 0.85, /ir = 0.05, £{T 0) is about equal to the 
largest system size at the lowest temperature, N — 2048. 
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